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ABSTRACT 

Context. Magnesium is an element of significant astrophysical importance, often traced in late-type stars using lines of 
neutral magnesium, which is expected to be subject to departures from local thermodynamic equilibrium (LTE). The 
importance of Mg, together with the unique range of spectral features in late-type stars probing different parts of the 
atom, as well as its relative simplicity from an atomic physics point of view, makes it a prime target and test bed for 
detailed ab initio non-LTE modelling in stellar atmospheres. Previous non-LTE modelling of spectral line formation has, 
however, been subject to uncertainties due to lack of accurate data for inelastic collisions with electrons and hydrogen 
atoms. 

Aims. In this paper we build and test a Mg model atom for spectral line formation in late-type stars with new or recent 
inelastic collision data and no associated free parameters. We aim to reduce these uncertainties and thereby improve 
the accuracy of Mg non-LTE modelling in late-type stars. 

Methods. For the low-lying states of Mgl, electron collision data were calculated using the A-matrix method. Hydrogen 
collision data, including charge transfer processes, were taken from recent calculations by some of us. Calculations for 
collisional broadening by neutral hydrogen were also performed where data were missing. These calculations, together 
with data from the literature, were used to build a model atom. This model was then employed in the context of 
standard non-LTE modelling in ID (including average 3D) model atmospheres in a small set of stellar atmosphere 
models. First, the modelling was tested by comparisons with observed spectra of benchmark stars with well-known 
parameters. Second, the spectral line behaviour and uncertainties were explored by extensive experiments in which sets 
of collisional data were changed or removed. 

Results. The modelled spectra agree well with observed spectra from benchmark stars, showing much better agreement 
with line profile shapes than with LTE modelling. The line-to-line scatter in the derived abundances shows some 
improvements compared to LTE (where the cores of strong lines must often be ignored), particularly when coupled with 
averaged 3D models. The observed Mg emission features at 7 and 12 pm in the spectra of the Sun and Arcturus, which 
are sensitive to the collision data, are reasonably well reproduced. Charge transfer with H is generally important as a 
thermalising mechanism in dwarfs, but less so in giants. Excitation due to collisions with H is found to be quite important 
in both giants and dwarfs. The A-matrix calculations for electron collisions also lead to significant differences compared 
to when approximate formulas are employed. The modelling predicts non-LTE abundance corrections AA(Mg) NLTE _ LTE 
in dwarfs, both solar metallicity and metal-poor, to be very small (of order 0.01 dex), even smaller than found in previous 
studies. In giants, corrections vary greatly between lines, but can be as large as 0.4 dex. 

Conclusions. Our results emphasise the need for accurate data of Mg collisions with both electrons and H atoms for 
precise non-LTE predictions of stellar spectra, but demonstrate that such data can be calculated and that ab initio 
non-LTE modelling without resort to free parameters is possible. In contrast to Li and Na, where only the introduction 
of charge transfer processes has led to differences with respect to earlier non-LTE modelling, the more complex case of 
Mg finds changes due to improvements in the data for collisional excitation by electrons and hydrogen atoms, as well 
as due to the charge transfer processes. Grids of departure coefficients and abundance corrections for a range of stellar 
parameters are planned for a forthcoming paper. 
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1. Introduction 

Neutral magnesium creates a broad range of spectral fea¬ 
tures in late-type stars, including some of the strongest lines 
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in their spectra. Consequently, Mg is detectable even in 
low-quality spectra and in nretal-poor stars, making Mg 
an excellent tracer of a-element abundances. The strong 
Mgl b lines with pressure-broadened wings can be used 
as a surface gravity diagnostic (e.g. lEdvardssonl ll988t 
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iFuhrmann et al.lll997t) . Additionally, the inter-combination 
line at 4571 A is unique among lines of similar strength 
in the Sun in that it is expected to have an almost LTE 
line_ formation across the entire profile (bee ICarlsson et ahl 
Il992l . and discussion therein). Emission features due to Mg I 
Rydberg transitions have been identified in the infrared 
(IR) reg ion of the solar spec t rum, particularly thos e near 
12 pm (jMurcrav et al.l [l98lb ICharig fc Novej fl 98. 'll) . The 
Mg I emission features dominate other features arising from 
similar transitions in other atoms due to a combination of 
abundance and ionisation energy. These emission lines have 
also been observe d in o ther m ain-seq uence and gian t sta rs: 
Arctu rus (lUitenbroek fc Noveslll996l), Pro cyon dRvde et all 
l2004 i). and Pollux ( Sundovist etT alIl2008fP1 . 

At T > 5000 K, Mg 1 is a minority species and as a 
result is expected to be sensitive to departures from LTE 
in stellar atmospheres, in particular due to photoionisation 
by non-local radiation. While in the Sun non-LTE effects 
on lines in the optical r ange are expected to be relatively 
weak (jZhao et ali 1199811 . significant departures from LTE 
have been predicted in metal-poor dwarfs and giants due to 
the increase in UV radiation becau se of the decreased line¬ 
blanketing dZhao fc Gehrenl l2000t) . Moreover, the 12 pm 
features have been satisfactorily explain ed by non-LTE 
modelling with a p hotospheric origin dCarlsson et ahlIl992b 
IZhao et. all 119981 : ISundovist et al.l I2fl08h . However! past 
studies have often been forced to use atomic collision data of 
questionable quality, frequently approximate formulae, and 
this is a si gnifica nt _source of uncertainty in the non-LTE 
modelling. ICarlsson et al.l showed that these uncertainties 
were the largest in their calculations of the solar 12 pm lines 
(see Sect. 6.3 of their paper). iZhao et ahl further demon¬ 
strated this, as well as showing the particular sensitivity 
of the 8806 A line to the rates for collisions with neutral 
hydrogen. 

The astrophysical importance of Mg, together with the 
unique range of spectral features in late-type stars probing 
different parts of the atom, plus its relative simplicity from 
an atomic physics point of view, makes it a prime target 
and test bed for detailed ab initio non-LTE modelling in 
stellar atmospheres. In this paper, we describe calculations 
of non-LTE Mg line formation based on a new model of 
the Mg atom with significant improvements in the collision 
data for neutral Mg. Our focus here is on inelastic colli¬ 
sion data, although some new data for collisional broad¬ 
ening are also presented. We perform calculations for ex¬ 
citation of the lower-lying levels due to electron impacts 
using the R-matrix method. Recent data for excitation 
and charge transfer due to hydrogen atom im pacts involv¬ 
ing low-lying levels calcu l ated by some of us (Guitou et al.1 


l201.lt iBelvaev et al.ll2012b iBarklem et al.ll2012 ) are also em¬ 

ployed. Furthermore, we have made considerable efforts 
to use physically motivated methods for calculating radia¬ 
tive and collisional data involving high-lying and Rydberg 
states. 

Additionally, in expectation of accurate stellar positions 
and motions from the GAIA satellite, several complemen¬ 
tary ground-based spectroscopic surveys have been initi¬ 
ated. Their goal is to determine high-precision abundances 
for eventually hundreds of thousands of stars and thereby 


1 The 12.3 pm feature is observed in absorption in Betelgeuse 
llUitenbroek fe Novedl 19961). b ut t his is du e to the blending with 
a strong water line llSundavist et alJl2008 l. 


unravel the formation and evolution of the Milky Way. The 
ongoing Gaia-ESO Survey will determine [Mg/Fe] ratios in 
all Galactic components to better than ~ 0.1 dex, in or¬ 
der to distinguish between the thin and thick disk and be- 
tween accreted and i n-situ halo com ponents, for instance 
(|Recio-Blancol [20131 : iGilmorel l2012f) . The targeted wave¬ 
length regions of the VLT/GIRAFFE spectrograph contain 
two strong Mg lines, 5528A and 8806A, detectable in all 
FGK type stars down to low metallicities. For stars in the 
solar neighbourhood, the GAL AH survey dDe Silva et al.1 
I2015T) aims at a higher precision (~ 0.05 dex) for many 
elem ents with the aim to identify dispersed star clus¬ 
ters (jZucker et al.l l2012f ). The wavelength settings of the 
AAT/HERMES instrument contain two weak Mg lines, 
571lA and 769lA, suitable for analysis of high-metallicity 
stars. 

The paper is structured as follows. In Sect. [2] we review 
previous studies and particularly the collision data used. In 
Sect. [3] we describe the model atom and present new data 
for inelastic electron collisions and collisional broadening 
due to hydrogen; in Sect. 14.11 the modelling is compared 
to observed spectra in benchmark stars. In Sect. 14.21 we 
compare our calculated profiles of the solar Mg I IR emission 
lines with observations. In Sect. Owe examine the effect of 
the different collisional processes in commonly used Mg 
lines. Finally in Sect. Owe present our conclusions. 


2. Previous works 

As mentioned, a significant number of non-LTE Mg I line 
formation studies have been conducted, and the main differ¬ 
ence between these, especially the more recent and compre¬ 
hensive studies since around 1990, and the work presented 
here lies with the atomic collision data for bound-bound 
processes. Thus, as a basis for comparison, in this section 
we present an overview of some past studies, focussing on 
the atomic collisional data employed. 

The earliest i n vestig ation of Mg I line formation by 
lAthav fc Ganfieldl (1969) examined the Mg b lines in the 
Sun. They used a small model atom with constant val¬ 
ues for the electron collisional cross sections based o n the 
van Regemorter (vR) formula (Ivan Regemorteilll962 ') - an 
often-used semi-empirical interpolation formula base d on 
the Bethe approximation. lAltrock fc Ganfieldl (1974) up- 
dated the co ll isiona l cross-sections to tho se calculated by 
Ivan Blerkoml (|1970H . iMauas et al.l (Il988fl used a 13-level 
model atom to study the 4571 A line and the Mg b line 
at 5173 A in the quiet Sun employing updated electron col¬ 
lisional data collected from various experiments and calcu¬ 
lations. 


Subsequent studies in the late 1980s and early 1990s fo¬ 
cussed on the Mg I IR emission features at 12 and 18 um 


observed in the Sun 

Lemke & Holwegeil 19871; [Chang et al. 

1991b ICarlsson et al. 

19921 among others), where high-lying 


and Rydberg levels were included. IChang et al.l used a 
model atom with 41 levels where 3s nl levels with l = s,p,d 
and n < 7 are included, and 7/, 7 q, 7 h, 7i are merged; 
8 < n < 15 are super levels. iGarlsson et al.l presented a 
model atom with 71 levels plus continuum, including up 
to 3s nl with n < 9, l < 8, while for n = 10, l = s,p 
levels are included separately and l > 2 are included as 
single super levels. For l > 3 singlet and triplet terms are 
collapsed to single levels. These studies demonstrated the 


Article number, page 2 of 1221 




















































































Y. Osorio et al.: Mg line formation in late-type stellar atmospheres: I. The model atom 


importance of comprehensive models including levels up to 
the quasi-continuum, where levels are collisionally domi¬ 
nated and strongly coupled to the continuum. This per¬ 
mits a population flow from the continuum to lower-lying 
levels. The importance of this was demonstrated by the 
failure of Lemke & Holweger’s model to reproduce the 12 
pm emission, due to a lack of levels sufficiently close to 
the continuum to be collisionally dominated, in addition 
to incorrect collision ra tes (see Sect. 6.6 of lCarlsson et all ). 
ISundovist et al.l (120081) also showed the need to extend the 


model atoms to include n > 10 for giant stars; the sparser 
atmosphere means that the levels become collisionally dom¬ 
inated at higher n. 

IChang et al.1 and ICarlsson et al. both too k electron col- 
lisional excitation rates from Mauas et all For the op¬ 


tically allow ed transitions not included in Mauas et al 


IChang et al.1 used the vR formula, while Ca rlsson et al.1 
used the impact param eter (IP) met hod (Seaton 19621) . 
For forbidden transitions, Cha ng et al.l used the vR formula 
with an arbitrary oscillator strength, while ICarlsson et al.l 
used an arbitrary scaling of nearby allowed transitions. 

Inelasticcollisions with neutral hydrogen atoms were in¬ 
cluded in IChang et al.l (ll99ll). For Mgl+ H collisional exci¬ 


tation the formula from lKaulakvsl (1198fil l was used for tran¬ 
sitions between Rydberg levels, but did not include data 
for the lower-lying Mg I levels. Charge-transfer processes 
involving Rydberg states and protons, Mg**(n/) + H + 
Mg + + H(n), are discussed, but suggested to be unimpor¬ 
tant in the solar photosphere, despite the large rate co¬ 
efficients, due to low number densities. We note that we 
have been unable to find the source of these data quoted 
as being in preparation. ICarlsson et al.1 discussed the pos¬ 
sible role of collisional excitation with neutral H, argu¬ 
ing that they will lead to significant 1-changing collisions 
among Rydberg states. They included this effect by set¬ 
ting high electron collision rates between levels of equal n 
and with small differences in quantum defects to ensure 
relative LTE be tween t hese l evels. Various revisions of the 
model atom in ICarlsson et al.l have been used for study¬ 
ing t he 12 and 18 pm features in stars other than the Sun 
(e.g. iRvde et ahl 1 2004 ISundovist et al.l120081) . some includ¬ 
ing collisions with hydrogen at oms v ia the Drawin formula 
(|Drawinl 119691 : [Steenbock fc Holwegeil 1 1 9841) . The Drawin 
formula has been shown in recent years to be unable to 
provide reliable estimates of hydrogen collision processes, 
and no estimates at all for optically forbi dden tr ansi tions 
or charge transf er processes (e.g. | Barklem et aI1l201ll) . 

The study of Zhao_et_a2j 1998h used similar energy lev¬ 
els to those of ICarlsson et al.1 some small differences and 
uncollapsed terms led to a 83-level atom. They adopted 
electron collisional data from the vR formula for all allowed 
transitions, while for forbidden transitions they assumed a 
collision strength of unity. Collisions with neutral hydrogen 
were included using the Drawin formula with an empiri¬ 
cally determined scaling factor that changes exponentially 
with the excitation energy of the upper level. Later studies 
by th e same group (|Gehren et al.l 12004: Ma shonkina et ahl 
l2008f) abandoned this exponential scaling in favour of a con- 
sta nt fact or. 

iPrzvbilla et al.1 (|2001l ) built an ato mic model that in¬ 
cludes the same levels as that of lCarlsson et al.l and added 
n = 11 and 12 for Mg I together with the Mgn levels 2 p 6 nl 
with n < 10, l < 4 and the ground state of Mg ill. Electron 
collisional excitation data for Mg I were calculated using the 


vR formula except for the transitions between the ground 
and the four lowest excited levels, which were taken from 
IClark et~ahl (119911) . No hydrogen collisions were included by 
them since, as they mentioned, hydrogen collisions will be 
unimportant in A-star s, the subject of that study. 

iMerle et~ahl (|201ll) used an updated model atom to 
study the non-LTE effects on lines to be observed by the 
GAIA satellite. The model atom included levels with nl 
quantum numbers up to n=10 and 1=9, including fine struc¬ 
ture at low l, leading to a 149 level model atom. They 
also included ~300 radiative inter-combination transitions. 
Electron collisional data were stated to be taken from 
“quantum mechanical” calculations, although the sources 
are not specified; the IP formula was used when calculations 
were unavailable. Hydrogen collisions were not included. 

Recently, iMashonkinal (|2013f ) performed tests on eight 
stars including the_ Sun u s ing a model atom with the 
same levels as I Zhao et~aLl (Il998i) for Mg I, 2 levels for 
Mg II and the ground state of Mg ill. The electron col¬ 
lisional data were the same as those of IZhao et ahl . but 
they replaced the rates for electron collisions from vR 
with th ose of Mauas et al.l when the latter were avail- 
able. IMa shonki nal us ed the hydrogen collision data from 
iBarklem et al. ( 20121) . 

Finally, we not e that all these studi e s used plane-parallel 
ID atmospheres. iRutten & Carlsson] (1 1994 s ! have investi¬ 
gated the effects of inhomogeneity in the 12 /im lines and 
demonstrated that the shapes of these lines are relatively 
insensitive to granulation-induced asymmetries. That work 
is, however, to our knowledge the only attempt at 3D 
non-LTE Mg line formati on to date. We further note that 
IRutten fc Carlssonl (|l994l ) provided an excelle n t sum mary 
of non-LTE mechanisms in Mg I, and iGhan 1 (11994 pro¬ 
vided a summary of the corresponding atomic physics re¬ 
lating to the infrared lines. 


3. Non-LTE modelling 

We performed standard ID non-LTE modelling, where 
the coupled radiative transfer and statistical equilibrium 
equations based on our Mg model atom are solved in 
ID stellar model at mospheres u sing version 2.3 of the 
MULTI code (ICarlssonlll986l Il992l). We employed the same 
background line opacity data (|Gollet et all 1 200.4 used in 
the calc ulation of the MARCS theoretical stellar atmo¬ 
spheres (Gustafss on et al.ll2008h . but re-sampled to 10 300 
frequency points. Turbulence was modelled by the usual 
micro- and macro-turbulence parameter^. 

The MULTI code treats the restricted non-LTE radia¬ 
tive transfer problem, in which the element of interest, here 
Mg, is treated as a trace element. Thus, it is assumed that 
there is no feedback on the atmospheric structure or back¬ 
ground opacities, which are computed assuming LTE. Only 
very few stud ies have treated the full non-L T E problem 
i n coo l stars (jShort fc Hauschildtl 120051 l2009h . iLind et ahl 
(|2012h considered the effects of increased electron densities 
and decreased UV opacity due to the increased ionisation 
of Fe I on the non-LTE line formation of Fe, arguing that 
the the first effect was insignificant and the second weaker 
than typical non-LTE corrections. However, such tests are 
limited, and as stated by the authors, only fully consistent 

2 Macroturbulence is assumed to be Gaussian and treated ac¬ 
cording to the radial-tangential model. 
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calculations will resolve this definitively. Such calculations 
are beyond the authors’ present capabilities and beyond 
the scope of the present work. We note that the abundance 
of Mg is almost always larger than Fe, especially in metal- 
poor stars (typically about a factor of 3), and Mg also has a 
somewhat lower ionisation energy, making it a more impor¬ 
tant electron donor. On the other hand, non-LTE does not 
perturb the ionisation balance of Mg as much as for Fe. It 
would be important to test the trace element assumption. 

In this section we present the details of our atomic 
model, often referred to as the “model atom”. In Sect cm we 
discuss the energy levels and the size of the model atom. 
In Sect. 1X^1 the radiative data are outlined, namely bound- 
bound (b-b) and bound-free (b-f) transition probabilities 
and line broadening cross sections, and in Sect. 13.31 the in¬ 
elastic collision data are described. In Sect. 13.41 we describe 
the modified model atoms used to test the influence of new 
data on the final results. 

3.1. Energy levels 

Following earlier studies, we decided to build a model atom 
including Mg I, Mg II, and the ground state of Mg ill. To 
build the structure of the model atom, e ner gy leve l d ata 
were collected from the NIST database (|Ralchenko et al.l 
1201 (t) . supplemented with data for high n (up to n = 80) 
in the case of Mg I from iKurucz & Belli (1 199511 . For high 
ang ular momentum levels absent from either the NIST or 
IKurucz & Belli collections, we made use of the Sommerfeld- 
Dirac expression to estimate the energies. In principle, it 
is desirable to have the model atom as large and as de¬ 
tailed as possible; however, this is neither computationally 
practical nor required to achieve converged results for the 
spectral lines of interest. Thus, to investigate the conver¬ 
gence of the model atom, a computer script was written to 
generate model atoms characterised by three parameters: 1) 
the largest principle quantum number for which levels are 
included, 2) the principle quantum number above which all 
levels are merged into “super levels” encompassing all states 
of same n, and 3) the term above which all fine-structure 
components J are merged into a single level. This permitted 
to test a range of atomic models in size and detail. 

Finally, after extensive testing, we found that the follow¬ 
ing model atom provided results that converged to better 
than 1% for a range of late-type stellar parameters. Mg I in¬ 
cludes states up to n = 20 and uses super levels for n > 10. 
Fine-structure splitting was implemented for 3 P states up 
to 3s Ip 3 P{ 2 , 3 , 4 }, as well as for the 3p 23 P{ 0 ,i,2} states. 
Splitting was not included for states 3 s nl 3 L where L > 2. 
Mg II includes states up to n = 11 and uses super levels for 
n > 7. For states up to 4/ 2 F| 5/ , 2 7 , 2 i fine-structure split¬ 
ting was included. This gives a final model atom covering 
Mg I, Mg II, with 108 and 34 states respectively, plus the 
ground state of Mg III. Grotrian diagrams for Mg I and n 
are shown in Fig. |T] 

3.2. Radiative data 

Radiative b-b t ransition proba b ilities were taken mainly 
from NIST and IKurucz & Belli (ll995lb When data for a 
transition were in both sources, the data from NIST were 
adopted. For the Mg b triplet, oscillator strengths were 
taken from lAldenius et~ah ( 2007 1. When transitions involv¬ 


ing Rydberg states of Mg I were needed and data were un- 
availa ble, tr ansition probabilities were calculated following 
iGuseinov &; Mamedovl (I2012H . The data for the most im¬ 
portant diagnostic lines are provided in Table [T| The final 
model atom has 1185 b-b transitions (962 Mg I and 223 
Mgll). In all cases, radiative damping was calculated con¬ 
sistently from the employed transition probabilities. 

For most UV lines, >B„ in the line formation region 
contributing to UV overionisation that tends to depopu¬ 
late lower levels. Thus, b-f processes are of considerable im¬ 
portance. Bound-free cross-sections for Mg I and Mg II were 
taken from TOPbase dGunto fc Mendozalll992tj Guiit'o et al. 
1993]). The data for Mg I were adopted from iButler et al~ 
(1199311 . and are available for n < 9 and l < 4. For Mgn 
the original source is unpublished and data are available 
for n < 10 and l < 3. For larger n and l , hydrogenic cross- 
sections were used. The threshold energies of the photoion¬ 
isation cross-sections taken from TOPbase do not coincide 
exactly with the experimental ones used in our modelling, 
and so the TOPbase cross-sections were shifted appropri¬ 
ately. Strong resonances in the b-f cross sections sometimes 
play an important role. The photoionisation cross-sections 
for the ground and first excited state of Mg I have no reso¬ 
nances strong enough to be important at typical tempera¬ 
tures in cool stellar atmospheres. The photoionisation from 
the 3s3p 1 P° level shows a resonance at A ~ 2900 A due 
to the 3 p 2 3 S auto-ionising level that is five times greater 
than the value of the cross-section at threshold, thus making 
it an important contributor. The photoionisation from the 
3s4s 1,3 S levels show resonances due to the 3p4s 1,3 P° auto- 
ionising levels that are more than 500 times greater than 
the value at threshold, making these resonances the main 
contributors to the total photoionisation rate at tempera¬ 
tures around 5000 K. Mgn photoionisation cross-sections 
do not have resonances in the energy range of interest. 

3.2.1. Line broadening 

In addition to turbulent and radiative broadening, both 
mentioned above, data on the effects of collisions on the 
modelled spectral lines were also required, in particular 
broadening of the lines due to collisions with charged par¬ 
ticles and neutral hydrogen atoms. For lines of diagnostic 
interest in this work (see table [T]), these data need to be 
as accurate as possible, while for other transitions this is 
probably not of great importance, and reasonable estimates 
suffice. 

For some of the Mg I lines we included collisional broad¬ 
ening due to electrons and protons using qu a dratic Stark 
broadening parameters from iDimitriievic fc Sahal-Brechotl 
(I1996T ). We adopted the line width at 5000 K ne¬ 
glecting variation with temperature, which is typi¬ 
cally w eak. For t he 18 | im l ines we adopted widths 
from Ivan Regemorter fc Hoang Binhl (Il993l) . For the inter¬ 
combination line 3 jP 3 P° — 3s 3 S at 4571 A we adopted the 
value used in iMashonkinal (f201 3f Fl. For all other lines (i.e. 
those not in table [TJ or of Mg Ii), Stark broadening was 
neglected. 

For transitions of Mgi involving low-lying states, col¬ 
lisional broadening due to collisions with neutral H is de¬ 
scribed via cross-sections and velocity parameters interpo- 

3 This value is stated to be taken from the VALD database, bit 
we were unable to find any such data for this line in the VALD. 
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Fig. 1: Grotrian diagrams of the Mg levels and transitions representing the model atom adopted for this study. Energy is 
relative to the ground state in each case. The nl configuration is shown for some levels. The dashed line in the diagrams is 
the ionised state. The vertical line in the 3 P term connects the 3s 3p 3 P° to the 3p 2 3 P level. Parity is not distinguished, 
and so these two terms are shown in the same column. 


lated in tables from the ABO theory (lAnstee fe O’Maral 
ll995tlBarklem k. Q’Maralll997tlBarklem et al.11998bU af). In 


the specific cases of the 4167 and 8806 A lines, these lie 
slightly outside the range of the tables and so specific cal¬ 
culations ha ve been made; the calculat ion for 8806 A was 
presented in iBarklem &; O’Maral (119971) . As such, calcula¬ 
tions are at the edge of the validity of the theory, and so 
these data must be expected to be more uncertain than 
for other lines. For transitions involving Rydberg states, 
the ABO theory is inappropriate for various reasons, most 
importantly becau se inelastic pr ocesses p la y a v ery impor¬ 
tant role. [Hoang Binh fc van Rcgcmorten (|1995l) . hereafter 
HBvR, have performed calculations for the 12 /tin lines 
in the impulse approximation. In this work we investigate 
these and several other lines in the infrared that involve 
Rydberg states (at 7 and 18 /im), and thus we performed 
broadening calculations for these lines using a similar for¬ 
malism described below. We show below that our results 
agree reasonably well with HBvR for the 12 /mi lines. 


In the impact approximation, the broadening cross- 
section for an isolated line correspondin g to the tran¬ 
sition i —> f can be written (e.g. iBarangerl ll962l : 
iHoang Binh fe van Regemorteiill995f) 



(1) 


If^n) - f f (n)\ 2 cm 


where and are inelastic cross-sections, and /) 

and ff are elastic scattering amplitudes. The broadening of 
the initial and final states due to inelastic scattering adds 
incoherently, and thus is simply the sum of the two contri¬ 
butions. The elastic scattering involving the two levels, how¬ 
ever, subtracts coherently and must be treated together at 
the scattering amplitude level. The inelastic cross-sections 


may for convenience be split into /-changing and n-changing 
collisions, such that 


inel 

J nl 


n'V 


Gnl^rn'V — / J ^nl^nl' H - 
V 


EE 

n'^n l' 


Gnl—tn'V 


( 2 ) 


The cross-sections for levels in transitions of inter¬ 
est were calculated and are presented in Table [U The /- 
changing inelastic cross-sections were calculated using the 
metdiod for inelastic Rydberg-hydrogen collisions due to 
IKaulakvsl <|1991h described in Sect. 13.3.21 The n-changing 
inelastic cross-sections were calculated in the scattering 
leng th approximation using the analytic expressions of 
iLebedev fc Marchenkol (119871 Eq. 25) and IKaulakvsl (119851 
eqns 4-5), which give practically identical results, as pointed 
out by Lebedev and Marchenko. We could calculate the n- 
changing cross-sections in the same way as done for the 
/-changing; however, the computational time is significant, 
and as seen in Table [21 the n-cha nging contr ibutio n is very 
weak, as expected (e.g. IHoang Binh fe van Regemorterl 
Il995l l. HBvR showed that the ratio of the elastic to the 
inelastic component is expected to be of order An/n 3 , and 
thus we did not include the elastic contribution. 

The total broadening cross-sections are then calculated 
according to Eq. [T]and are presented in Tabled The cross- 
sections for a collision velocity of 10 4 m/s are presented, 
together with a velocity parameter a that was calculated 
by fitting a powe r-law behaviour of the c ross-section with 
velocity following lAnstee fc O’Maral (Il995fl . Line widths are 
also presented and compared with the results of HBvR for 
the 12 fi m lines. As seen from Tabled our cross-sections 
are generally larger than those of HBvR, although notciuite 
as large as their calculations using the method of lKaulakvsl 
(jl99lf ) . Note that we have included the n-changing com¬ 
ponent, which was neglected by HBvR. Our cross-sections 
for the 12.2 and 12.3 /in 1 lines differ from theirs by 1 and 
21 per cent, respectively. The line widths differ by consid¬ 
erably more because HBvR used an approximate method 
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Table 1: Data for important transitions, including diagnostic lines. Columns 1 and 2 give the transition wavelength and 
labels for the states involved. Columns 3 and 4 give the adopted oscillator strength, in the form log(g/), and the source 
of these data (see notes). Columns 5-8 give collisional broadening data (see text and notes). Columns 9-12 show, for each 
line of Mg I, the collision process most responsible for the changes between model B and model F, and thus indicates 
the most important contributor of the new collision data introduced in this work (see Sect. [5] for more details). Results 
are given for the four test atmospheric models described in Sect. [5j Dr = Dwarf, metal-rich; Dp = Dwarf, metal-poor; 
Gr = Giant, metal-rich; Gp = Giant, metal-poor. The collisional processes are labelled (see Sectj5]) H when the line has 
a sensitivity similar to charge transfer with H and CHO; collisional excitation with H, CH, and CE indicates the line is 
sensitive mostly to the electron collision rates. For Mg II lines there are non-LTE effects in some lines, but the non-LTE 
abundance corrections are not sensitive to any particular collisional process. 


A 

Transition 

log gf 

Source T 6 

log r 4 /N e Source 

Dr Dp 

Gr 

Gp 




£7 ot 

(5000 K) 




[Al 



[a.u.] 

[rad s —1 cm 3 ] 





Mg i 


3 829 

3d 3 D - 3p 3 Pq 

-0.23 

Fro/NIST 

708 

3 832 

3d 3 D -3p 3 P° 

0.25 

Fro/NIST 

708 

3 838 

3d 3 D -3p 3 P 2 

0.47 

Fro/NIST 

708 

4 167 

Id 1 D - 3p 1 P° 

-0.75 

C-T/NIST 

222 

4 571 

3p 3 Pi — 3s 1 S 

-5.62 

Fro/NIST 

222 

4 703 

3d 1 D - 3p 1 P° 

-0.44 

C-T/NIST 

2806 

5 167 

4s 3 S - 3p 3 Pq 

-0.93 

Aid 

728 

5 173 

4s 3 S-3p 3 P° 

-0.45 

Aid 

728 

5 184 

4s 3 S-3p 3 P° 

-0.24 

Aid 

728 

5 528 

4d 1 D-3p 1 P° 

-0.50 

C-T/NIST 

1460 

8 710 

7d 3 D - 4p 3 Pq 

-1.57 

But/NIST 


8 713 

7d 3 D-4p 3 P° 

-1.09 

But/NIST 


8 718 

7d 3 D-4p 3 P 2 

-0.87 

But/NIST 


8 736 

7/ 3 F° — 3d 3 D 

-0.53 

But/NIST 


8 806 

3d 1 D-3p 1 P“ 

-0.13 

Fro/NIST 

529 

73 700 

6 h H°-5g 3 G 

1.34 

Civ 

4950 

122 200 

7 h H° — 6 p G 

1.62 

Civ 

5191 

123 200 

7i I - 6 /r H° 

1.95 

Hydro 

4657 

188 300 

8 h H° — 7p G 

0.49 

Hydro 

4497 

189 500 

Si l -7h H° 

1.90 

Hydro 

4304 


Mg ii 




4 385 

bd 2 D - 4 p 2 P °/2 

-0.78 

Sie/NIST 


4 391 

5d 2 D-4p 2 P ° /2 

-0.48 

Sie/NIST 


4 481 

4/ 2 F° - 3d 2 D 

0.76 

Fro/NIST 


5 402 

7p 2 G-4/ 2 F° 

0.06 

K-P 


7 877 

4d 2 D 3/2 -4p 2 P ° /2 

0.39 

Sie/NIST 


7 896 

4d 2 D 5/2 -4p 2 P ° /2 

0.64 

Sie/NIST 



0.301 

-4.51 

D-SB 

H 

H 

CH 

CH 

0.301 

-4.51 

D-SB 

H 

H 

CH 

CH 

0.301 

-4.51 

D-SB 

H 

H 

CH 

CH 

0.249 

-3.49 

D-SB 

H 

H 

CH 

H 

0.249 

-6.51 

Mas 

H 

H 

CE 

CE 

0.269* 

-4.11 

D-SB 

H 

H 

CH 

H 

0.238 

-5.40 

D-SB 

H 

CH 

CH 

CH 

0.238 

-5.40 

D-SB 

CH 

CH 

CH 

CH 

0.238 

-5.40 

D-SB 

CH 

CH 

CH 

CH 

0.312 

-4.56 

D-SB 

H 

H 

CH 

CH 


-2.72 

D-SB 

CHO 

CHO 

CH 

H 


-2.72 

D-SB 

CHO 

CHO 

CH 

H 


-2.72 

D-SB 

CHO 

CHO 

CH 

H 


-2.95 

D-SB 

CHO 

CHO 

CH 

H 

0.277 

-5.39 

D-SB 

CH 

CH 

CH 

CH 

1.549* 

-3.06 

vR-HB 

CH 

CE 

CH 

CH 

1.738* 

-2.39 

D-SB 

CH 

CH 

CH 

CH 

1.752* 

-2.55 

D-SB 

H 

CHO 

CH 

CHO 

1.764* 

-2.92 

vR-HB 

CHO 

CHO 

CH 

CH 

1.778* 

- 2.12 

vR-HB 

H 

H 

CH 

CH 


Notes. Oscillator stre ngths (/-values) were collected m ostly f rom the NIST databas e (Ralchcnko ct al. 201 (3). The ori ginal sources 


of the NIST data are lFroese Fischer fc Tachievl (2010) [Fro|, 
(1 19Q8 Ti [Siej. Other sources are [Aldenius et all lt2007li [Aid], 


Chang^JTand Tj_990|) [C-TJ Butler et al.l (1 1 9 031) [Butl, Siegel et all 
Civis. S. et al.l ll 201 .'ll [Civ] and iKurucz fc Pcvtremann ( 1973) [K- 


[Civ] and IKurucz fc Pevtremaim ~ lfl97^r [K- 
P[. Fo r the Mg l IR. line s wi th no data found in the literature we calculated the /-values using the hydrogenic formula from 
iGuseinov fc Mamedov! (1 2012 ) [Hydro]. To calculate the van der Waals line widths Tg we used data in the ABO theory format 
where o is the broadening cross-section in atomic units and a is the velocity parameter. These were calculated with the ABO 
theory except for those marked with an asterisks, which were ca l culat ed in this work (see Table[3]for the IR line s ). Sta rk broadening 
line widths, T 4 , we re taken from [Dimitriievic fc Sahal-Brcchot (1996) [D-SB], van Regemo rter fe Hoang Binbl (1993) [vR-HB] and 
llMashonkinall20i3i . including only electrons) [Mas], 


for calculating the line widths from a cross-section at a sin¬ 
gle velocity, rather than an integration over a Maxwellian 
velocity distribution. In any case, we conclude that the un¬ 
certainty in the broadening cross-sections for these lines 
involving Rydberg states is at least 20 per cent, probably 
larger due to approximations inherent to both methods. 


For transitions covered neither by the ABO theory nor 
by the calculations for infrared lines described above, the 
Unsold formula was used with an enhancement factor of 
2.5. This includes all transitions of Mg II. 


3.3. Inelastic collision data 

In this subsection we describe the data adopted for inelastic 
collision processes: collisional excitation and ionisation due 
to electrons, and collisional excitation and charge transfer 
processes due to hydrogen atoms. This includes the descrip¬ 
tion of new calculations for collisional excitation of Mgl 
using the /(-matrix method. 

Figure [5] visualises the total downward collision rates 
calculated in the (3Dq)id atmospheric model (see Sect. 14. l(l 
at a depth corresponding to roughly T = 6000 K 
(log , r 500nm =—0.23) using our final model atom. The collision 
matrix fo r Mg I has a structure very similar to that seen in 
ICarlsson et al.l (|1992[ ). In particular, we see the block struc¬ 
ture of the very strongest transitions along the diagonal 
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Table 2: Inelastic /-changing QT/, a) and n-changing (^ , a) cross-sections for levels of interest, in atomic units. Columns 
2 and 3 give results from this work at v = 10 4 m/s. For direct comparison with /-changing crosssections from Table 4 
of HBvR, results are also given for two states at v = 10 457 m/s corresponding to the average collision velocity at 
T = 5000 K. Column 4 gives our results and Col. 5 the results from HBvR. 


Level 

V = 

10 4 m/s 

V = 

10 457m/s 

Ei'0- 

[au] 

[au] 

b _ 

d 1 

(BBvR) 

[au] 

5g 

4752 

78 



6g 

5456 

130 



6h 

4938 

132 

4553 

4036 

7g 

4744 

196 



7h 

4600 

196 



7i 

4048 

196 

3714 

3088 

8h 

3786 

270 



8i 

3542 

270 



8k 

3146 

270 




Table 3: Broadening cross-sections, velocity parameters, and line widths for the lines of interest. Columns 3 and 4 give 
broadening parameters, the cross-section <r br at v = 10 4 m/s, and the velocity parameter a (the ABO theory format). 
Columns 5 and 6 compare the cross-sections from this work with those of BBvR at v = 10457 m/s. The broadening 
cross-section can be determined directly from their line width ((7 h/A ? h)/2 x 1/(wh) hi their notation). In Cols. 7 and 8 
the line widths (half-width at half-maximum) at T = 5000 K are compared. 


Wavelength 

Transition 


a 

cr br (10 457m/s) 

w/Nn 

w 


[au] 


[au] 

(BBvR) 

[au] 

[cm 3 

(BBvR) 
rad s 1 ] 

7 

5g — 6 h 

4950 

1.549 



14.9 


12.2 

6g - 7 h 

5191 

1.738 

4804 

4840 

16.3 

14.2 

12.3 

6 h — 7 i 

4657 

1.752 

4303 

3560 

14.7 

10.5 

18.8 

7g-8h 

4497 

1.764 



14.3 


18.95 

7 h - 8 i 

4304 

1.778 



13.7 


18.96 

7 i - 8 k 

3831 

1.789 



12.2 



caused by the strong /-changing collisions where An = 0, 
where collisions with Al = 1 are strongest of all. We also see 
structure off the diagonal due to the strength of collisions 
with Al = 1 and An = 1. A new feature in this work is 
the structure seen in the collisional coupling of Mg I states 
to the Mg II ground state. Inclusion of charge transfer pro¬ 
cesses leads to rates between the levels around 4s 4 S and 
the Mg II ground state having the same order of magnitude 
as the collisional rates from these levels to neighbouring 
levels. 

Before we describe the data, a word on the 3 p 2 3 P level is 
warranted. This is the only level included in our modelling 
where the dominant configuration is not 3s nl. This level 
is too highly excited to be included in the detailed quan¬ 
tum mechanical calculations, and it is questionable whether 
the approximate methods are applicable to transitions in¬ 
volving two electron^. Our expectation is that collisional 
processes involving this state will be relatively weak, and 
thus we ignored collisional coupling to this state (as can be 
seen in Fig. [5]). This causes this level to behave quite dif¬ 
ferently from other levels of Mg I in the non-LTE modelling 
and to depart from equilibrium very deep in the atmosphere 


4 The 3s3p 1,3 P° — 3p 2 3 p transition involves only one electron 
and is an exception. However, we still neglect the collisional 
coupling between these levels, which can be expected to be weak 
as a result of the large energy difference. 


(logT 500nm ~ 1 in the solar case). Bowever, this is not im¬ 
portant for any lines observed in stellar spectra. 

3.3.1. Electrons 

Collisional excitation and de-excitation between low-lying 
states of Mg I are accounted for with //-matrix calcula¬ 
tions, which are presented in detail below. These calcula¬ 
tions provide data for the transitions between the lowest 
ten states of Mg I, up to 5p 3 P°. For transitions involving 
higher lying states, including Rydberg levels, we must use 
more approximate methods. For allowed transitions, two 
such methods a re in common use in stellar astrophysics: 
the vR formula (|van Regemortei(ll962l f. which is based on 
the Born and Bethe approximations with an empirical ef¬ 
fective Gaunt factor, a nd t he s emi-cl assical impact param¬ 
eter (IP) method from ISeatonl (Il962h . The IP method has 
been extended to the case of positive ions, where electrons 
have hyperbolic tr aject or ies i n the semi-classical picture, 
bv iBurgess & Summerd ( 1976 1. but this seems to be little 
used. 

It is well known that the Born approximation overes¬ 
timates cross-sections at low energies and the IP method 
gives better resul ts, as was repor ted f or inst ance bv ISeator] 
(Il96l : see also iBelv fe van Regemorted (11970P for a dis¬ 
cussion. Figure El compares results from the two methods 
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Table 4: Energies for Mgi states from electronic structure 
calculations and comparison with experimental values. 


State 

Energy (eV) 
Experiment CIVPOL 

2p 6 3s 2 

3 S 

— 

— 

3s3p 

3po 

2.71 

2.67 

3s3p 

lp o 

4.34 

4.38 

3s4s 

3 S 

5.11 

5.10 

3s4s 

3 S 

5.39 

5.36 

3s3d 

X D 

5.75 

5.73 

3s4p 

3po 

5.93 

5.88 

3s3d 

3 d 

5.94 

6.01 

3s4p 

lp o 

6.12 

6.05 

3s5s 

3 S 

6.43 

6.43 


Notes. The energies shown are relative to the ground state 
of Mg l. The experim ental values were taken from NIST 
dRalchenko et ail l ^OlOT ). 


Fig. 2: Comparison of the Mgl electron collision rate co¬ 
efficients at T=5000 K using the impact parameter (IP) 
method and the van Regemorter formula (vR). 


for allowed transitions in Mgl. The vR results are larger 
than those from IP. We also performed tests using either 
IP or vR in the atomic model, and we found that while we 
can reproduce the IR Mgl emission lines observed in the 
Sun when using the IP method, we could not do so when 
using the vR formula (see Fig. [5]). Thus, there is both phys¬ 
ical and astrophysical support for preferring the IP method 
data over vR when i?-matrix calculations are not available. 

Collisional ionisation of low-lying states by electrons was 


implemented according to the hydrogenic approximation 
presented in lCoxI (120001 Sect 3.6.1), which is based o n sem i- 
empir ical expressions taken from Bely fe v an Regemorted 
(|1970H . which originate from lPer cival| (|1966). For R ydberg 
levels we used the formula from IVrinceanu (|2005fb where 
the dependency of the collisional cross-sections on the total 
electronic orbital angular momentum of the Rydberg atom 
is taken into account. 


i?-matrix calculations. We now give the details of the elec¬ 
tron collision calculations for the lowest ten states of Mg I 
using the R-matrix (RM) method. All calculations were 
performed in L-S coupling. Electronic structure_ cal cula¬ 
tions were pe rformed using the cod e CIVPOL (Hib bert et al.1 
see also lPlummer et alll20()4h . which employs the con¬ 
figuration inte racti on method. CIVPOL is an adaptation of 
CIV3 (iHibbertJl 1975 1. allowing the construction of polarisa¬ 
tion pseudo-states. Orbitals for core electrons, Is, 2s, 2 p 
and 3s, were taken from the gr o und-s tate Hartree-Fock 
calculations of lClementi fc R.oettil (119741 ). Valence orbitals, 
4s, 5s, 3 p, 4p, 3d, 4 d and 4/, were optimised with CIVPOL, 
providing energy levels that agree well with experimental 
values taken from NIST; see Table [|] 

We calculated the excitation cross-sections with the 


RM method (B urke et, al.1 Il97lb iBurke fc Robbl 119761 : 
iBerrington et all 1974 . 197811 . This method separates the 
calculation into two cases: when the colliding electron is in 
the internal region, that is, close to the target atom where 


interactions are strong, and when the colliding electron is in 
the external region, far enough from the target atom such 
that interactions are weak and correlation and exchange 
effects can be neglected. The A-matrix ensures continu¬ 
ity between wave functions at the boundary between the 
two regions. To perform the RM calculati ons fo r t he inter¬ 
nal r egion we used the code RMATRIX I (IBerrington et all 
119951) , and for the external region a version of the code 
STGF (IBerrington et all 1 1987T) modified to treat collisions 


i 


with neutral atoms (lBadnelllll999l f. 

After exploring the effects of the input data on the final 
collisional cross-sections, we adopted a target atom consist¬ 
ing of 25 spectroscopic states (i.e., energy levels), with ten 
additional pseudo-states states to represent the Rydberg 
levels and the continuum, giving 35 states in total. Target 
atoms with more levels were also tested with no significant 
changes in the resulting cross-sections. 

The maximum occupation numbers in each nl orbital 
used to build the configurations for both N and N + 1 
atomic structures were 


Orbital 
N electrons 
N + 1 electrons 


3s 4s 5s 3 p 4p 3d 4 d 4/ 
22121211 
2223232 2. 


The boundary between the internal and external regions 
was set at 30.4 atomic units. The final calculation included 
partial waves L < 20. We made test calculations using par¬ 
tial waves up to L < 30 to check convergence. i?-matrix 
calculations using L < 10 provided rate coefficients within 
20% of those calculated using L < 20 except for the 3d 3 D - 
4p 3 P° transition, where the rate coefficient increased by a 
factor of 2. There were no significant differences between the 
resulting rate coefficients when using L < 30 and L < 20 in 
the temperature range typical of cool stellar atmospheres. 

The final collisional cross-sections obtained were com¬ 
pared wit h experiments an d with other calculations where 
available. IZatsarinnv et al.l ( 2009 ) calculated electron col¬ 
lisional excitation cross-sections for Mgl from the ground 
state 3 s 1 S to the excited states 3p 1,3 P°, 3d 1 D, 4s 3 S and 
4 p ^R 0 and compared them with other calculations and 
with experiments. The cross-sections calculated in this work 
agree with those within ~20%. Rate coefficients were ob¬ 
tained by integrating the final cross-sections weighted by a 
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Fig. 3: Comparison between the electron collisional rates 
at T = 5000 K calculated by us and those obtained with 
the impact parameter method (IP, diamonds) and the van 
Regermorter formula (vR, triangles). 


Fig. 4: Comparison between the electron collisional rates at 
T = 5 000 K calculated by us and those used bv lMauas et all 
(|1988ll . Diamonds corresp ond to the a llowed transitions 
used by them (Table 1 in iMauas et al.ll an d trian gles are 


the forbidden transitions (Table 2 in IMauas et ali i. 


Maxwellian velocity distribution and are presented in Ta¬ 
ble [Gl 

Figure [3] compares our results with those obtained with 
the IP and vR formulas. Both approximate formulas typ¬ 
ically yield lower values than the RM results. Most past 
studies of Mg in non-LTE have, however, employed the elec¬ 
tron collisional data for transitions be t ween l o w-lying levels 
that were presented in IMauas et al.1 (119881) . IMauas et al.l 
studied Mg I non-LTE effects in the Sun using a 13-level 
model atom, and their compilation of electron collision 
data, adopting data from different sources (see that pa¬ 
per for details), has been employed extensively since. Thus, 
it is of interest to compare our results with their compila¬ 
tion, which is done in Fig. 0] Allowed transitions (circles) 
have results similar to ours, but for forbidden transitions 
(diamonds) our res ults are s ignificantly larger. For forbid¬ 
den transitions IMauas et al.i mostly used rates equal to 0.1 
times the collisional rate of an allowed transition with a 
similar energy difference calculated using the vR formula. 

For optically forbidden transitions, there are no sim¬ 
ple approximate formula similar to vR or IP, and we re¬ 
quire some method to estimate the rates of forbidden tran¬ 
sitions between highly excited levels. Common approaches 
to circumvent this problem are to adopt an arbitrary col¬ 
lision strength fly of the order of unity, or to assume col¬ 
lision strengths similar to nearby optically allowed transi¬ 
tions. We adopted a more nuanced approach. First, from 
a physical point of view, it makes sense to divide opti¬ 
cally forbidden transitions into two groups: those involving 
electron exchange (i.e. corresponding to intersystem tran¬ 
sitions with change of spin state) and those not involving 
electron exchange. At high energy, transitions not involv¬ 
ing electron exchange may be described in the Born ap¬ 
proximation, while transitions involving electron exchange 
are described in the Ochkur approximation. This leads to 
significantly different behaviour of the cross-sections with 


at least at high ene rgy (e.g. iBransden fc .Toachainl 
Burgess fc Tullv 199^ 1. Second, we then assume that 
transitions of the same type (exchange or non-exchange) 
and same transition en ergy will have similar r ate coeffi¬ 
cients. Using Eq. 20 of I Burgess & Tulivl (ll992lh this im¬ 
plies a constant value of T ij/gi at a given temperature, 
where T,;, is the therm ally averaged collision strength (e.g. 
Burgess fc Tulivl Il992. eqn. 21), and gi is the statistical 
weight of the initial state i. 


Indeed, we found using the results of the RM calcula¬ 
tions that the behaviour of Yy / gi for a given transition 
type was similar, but significantly different between the 
two groups. Among the calculated RM collisions we have 

4 non-exchange forbidden transitions and 40 spin-exchange 
transitions. We calculated the mean value for T ij/gi for 
the two groups from the RM calculations, and the results 
are presented in Table [5] The scatter of the spin-exchange 
transitions about the mean value is greatest at low tem¬ 
peratures, becoming smaller at higher temperatures. At 

5 000 K the values span over three orders of magnitude 
(0-01 Ty/g.j < 10) with half of the values within a factor 
of 2 of the mean. The mean behaviour for exchange transi¬ 
tions with temperature is rather flat (except at the lowest 
temperature), while for non-exchange transitions the values 
steadily increase over the temperature range of interest. To 
estimate rate coefficients for optically forbidden transitions 
not included in the RM calcu lations, we inter polate in this 
table, and then use Eq. 20 of lBureess fe Tulivl (Il992f) . which 
relates the thermally averaged collision strength to the rate 
coefficient. 
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Table 6: Rate coefficients for de-excitation of Mgi by electron collisions. Rate coefficients, ( av ), are given in cm 3 s 


3 0 —1 


\ final On lc 

initial \ ^ 

3p 3 P° 

3p i P° 

4s 3 S 

4s i S° 

3d i D 

4p 3 P° 

3d 3 D 

4p 1 P° 





1000 K 





3p 3 P° 

5.01E-09 









2, v l V° 

3.95E-08 

2.18E-06 








4s 3 S 

4.25E-09 

7.69E-07 

5.13E-08 







4s 3 S 

1.82E-09 

2.59E-07 

9.94E-08 

7.94E-08 






3d 1 !) 

5.52E-10 

8.21E-08 

2.58E-08 

6.35E-08 

1.25E-08 





4p 3 P° 

1.72E-09 

5.69E-07 

4.73E-08 

5.42E-07 

1.13E-07 

1.73E-06 




3d 3 D 

3.05E-10 

7.52E-08 

6.83E-09 

7.95E-08 

1.56E-08 

1.64E-07 

8.61E-05 



4p 1 P° 

1.58E-09 

1.24E-07 

4.14E-08 

5.18E-08 

1.64E-07 

3.29E-06 

1.23E-07 

4.94E-07 


5s 3 S 

1.56E-06 

3.11E—04 

2.56E-05 

2.51E-04 

1.17E-04 

1.19E-03 

9.01E-05 

1.11E-03 

1.55E-04 





2 000 K 





3p 3 P° 

1.79E-08 









3p 1 P° 

5.65E-08 

6.79E-07 








4s 3 S 

1.11E—08 

4.57E-07 

7.91E-08 







4s 3 S 

7.37E-09 

1.74E-07 

2.22E-07 

1.31E-07 






3d 3 D 

2.77E-09 

9.74E-08 

8.36E-08 

9.25E-08 

2.47E-08 





4p 3 P° 

3.10E-09 

2.68E-07 

4.95E-08 

3.98E-07 

5.87E-08 

6.24E-07 




3d 3 D 

8.93E-10 

6.68E-08 

1.53E-08 

9.98E-08 

1.45E-08 

1.22E-07 

4.25E-05 



4p 1 P° 

4.24E-09 

8.25E-08 

6.73E-08 

4.27E-08 

1.93E-07 

1.70E-06 

1.40E-07 

2.61E-07 


5s 3 S 

1.10E-07 

5.30E-06 

1.21E-06 

5.57E-06 

2.58E-06 

1.98E-05 

4.17E-06 

2.26E-05 

6.83E-06 





5 000 K 





3p 3 P° 

3.18E-08 










6.71E-08 

3.04E-07 








4s 3 S 

1.54E-08 

2.55E-07 

8.09E-08 







4s 3 S 

1.45E-08 

1.18E-07 

3.67E-07 

1.42E-07 






3 d 3 D 

6.50E-09 

8.58E-08 

1.26E-07 

7.86E-08 

3.85E-08 





4p 3 P° 

3.85E-09 

1.62E-07 

4.37E-08 

3.08E-07 

3.08E-08 

2.87E-07 




3d 3 D 

1.38E-09 

5.31E-08 

1.98E-08 

1.42E-07 

9.67E-09 

8.60E-08 

2.24E-05 



4p 1 P° 

8.14E-09 

6.72E-08 

9.58E-08 

3.48E-08 

3.10E-07 

1.51E-06 

1.33E-07 

1.82E-07 


5s 3 S 

2.02E-08 

4.21E-07 

1.99E-07 

4.20E-07 

2.57E-07 

1.84E-06 

6.18E-07 

1.81E-06 

1.36E-06 





8 000 K 





3p 3 P° 

3.19E-08 









?yp 1 P° 

6.93E-08 

2.25E-07 








4s 3 S 

1.36E-08 

1.95E-07 

6.84E-08 







4s 3 S 

1.51E-08 

9.31E-08 

4.21E-07 

1.24E-07 






3d 3 D 

7.81E-09 

7.51E-08 

1.25E-07 

5.98E-08 

4.14E-08 





4p 3 P° 

3.65E-09 

1.41E-07 

3.92E-08 

2.69E-07 

2.16E-08 

2.04E-07 




3d 3 D 

1.38E-09 

4.53E-08 

1.84E-08 

1.85E-07 

7.41E-09 

6.64E-08 

1.66E-05 



4p 1 P° 

1.02E-08 

6.05E-08 

1.12E-07 

2.98E-08 

4.13E-07 

1.55E-06 

1.13E-07 

1.56E-07 


5s 3 S 

1.33E-08 

2.25E-07 

1.52E-07 

1.83E-07 

1.54E-07 

1.03E-06 

3.47E-07 

8.21E-07 

1.03E-06 





10 000 K 





3p 3 P° 

3.04E-08 









3p 1 P° 

7.05E-08 

1.95E-07 








4s 3 S 

1.22E-08 

1.73E-07 

6.09E-08 







4s 3 S 

1.49E-08 

8.22E-08 

4.44E-07 

1.13E-07 






3d 3 D 

8.22E-09 

6.96E-08 

1.23E-07 

5.10E-08 

4.19E-08 





4p 3 P° 

3.46E-09 

1.35E-07 

3.68E-08 

2.51E-07 

1.81E-08 

1.74E-07 




3d 3 D 

1.33E-09 

4.21E-08 

1.72E-08 

2.10E-07 

6.50E-09 

5.75E-08 

1.44E-05 



4p 1 P° 

1.13E-08 

5.70E-08 

1.21E-07 

2.73E-08 

4.70E-07 

1.56E-06 

1.02E-07 

1.44E-07 


5s 3 S 

1.16E-08 

1.80E-07 

1.49E-07 

1.31E-07 

1.34E-07 

8.60E-07 

2.75E-07 

5.98E-07 

9.61E-07 


3.3.2. Hydrogen atoms 

The rate coefficients, (crv), for excitation and de-excitation 
processes, 

Mg(3s nl 2S+1 L) + H(ls) ^ Mg(3s n'l' 2S ' +1 L') + H(ls), 
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Table 5: Average values of T,y / g t (dimensionless) calculated 
from RM calculations as a function of temperature T. 


T |K| 

non-exchange 

exchange 

1000 

0.086 

1.020 

2 000 

0.174 

0.654 

3 000 

0.257 

0.635 

4 000 

0.336 

0.642 

5 000 

0.413 

0.650 

6 000 

0.492 

0.655 

7 000 

0.568 

0.658 

8 000 

0.646 

0.660 

10 000 

0.797 

0.658 


(3) 

and for the charge transfer processes, ion-pair production 
and mutual-neutralisation, involving the ionic state, 

Mg(3s nl 2S+1 L) + H(1s) ^ Mg+(3s 2 S) +H", (4) 

are taken fromjBajklgm_etal. _j [2012f) . which are based on 
calcul a tions of lBelvaev et al.l " 20121) : see also iGuitou et al.l 
(l20TT)l farnl) . These ab initio quantum chemical calcula¬ 
tions on which the scattering calc ulations are ba sed were 
obtained using the code MOLPRO (IWerner et al.l 120091) and 
cover the seven lowest lying states of Mg I and the ground 
state of Mg ii. 

We wish to emphasize that with regard to hydro gen 
collis i ons in general, the accurate quantum fBarklcm et al.l 
120031 : IBarklem et al.l 12010L 120121) as w ell as the more ap¬ 
proxi mate model calculations (|Belvaevll2013l : iBelvaev et al.l 
i2on . based on a correct physical background, provide non¬ 
zero rate coefficients for both optically allowed and op¬ 
tically forbidden transitions, as well as for charge trans¬ 
fer processes, which generally have the highest rate coef¬ 
ficients among all hydrogen collision processes. Thus, the 
Drawin formula fails to provide reliable estimates for im¬ 
portant inelastic processes in collisions with hydrogen. In 
th e case of Mg, the quantum eight-channel calculations 
IBelvaev et al.l ( 20121 ) provide 28 non-zero cross-sections for 
the endothermic processes (excitation and ion-pair produc¬ 
tion), while the Drawin formula gives only five non-zero 
(ove rest imat ed) cro ss-sec tions for the same 28 transitions 
(see IBarklem et ahll2012D . This is a general feature of the 
Drawin formula, as discussed in detail by Barklem et a l. 

mm- 

For collisional excitation involving higher lying and Ry¬ 
dberg _states of Mg I, we e mployed the free electron model 
of IKaulakvsl (11986 , 119911) , which considers a binary en¬ 
counter of the Rydberg electron with a neutral perturber. 
This approach allows writing the cross-section for the in¬ 
elastic transition nl —» n'V in terms of the elastic differen¬ 
tial cross-section for the electron-perturber scattering and 
th e m ome ntum- space wave functions of the initial state 
(IKaulakvslIl991l eqs. 6-8). For the wave functions g n i(p), 
where p is the electron momentum, we calculated non- 
hydro genic wave functions fr om q uantum-defect theory fol- 
lowing fHoang Binh fe van Regemorteil (|1997f Fl. For neutral 


5 The required Hankel transforms were derived analyti¬ 
cally for l ranging from 0 to 13 using Mathematica. 
We note that the denomina t ors in Eqns. 32 and 33 of 
iHoang Binh fe van Regemorteil (119971 ) contain a misprint where 
an exponent is incorrectly written as a multiplicative factor - 
the denominator of Eqn. 32 should be [u~ 2 + q p]("-*+ 1 )/ 2 . 


hydrogen perturbers, the relevant elastic differential cross- 
section \f e {p, 9)| 2 is that for e _ +H collisions, where f e (p, 9) 
is the scattering amplitude. This total cross-section can be 
written in terms of the singlet and triplet amplitudes, f + 
and /“, which is found by averaging over initial spins and 
summing over final spins to obtain the well-known expres¬ 
sion (e.g. ISeateinlil962j) 

\fe(p,9)\ 2 = i|/+(p,(9)| 2 + ^\f~(p,9)\ 2 . (5) 

The elastic singlet and triplet scattering amplitudes for 
e _ +H collisions are known from variational calculat ions of 
the phase shifts for s- and p -wave scattering (|Schwart,dl96ll : 
lArmsteadlll968l) . We make use of the form introduced by 
HBvR, namely 

\f ± (p,9)\ 2 =A ± +B ± cos(9), (6) 

where A ± is the s-wave contribution and the s-p cou¬ 
pling. The coefficients A± and B ± are fisted as a function 
of p in Table 3 of HBvR, derived from the published phase 
shifts mentioned above. From this we calculated the total 
cross-section for the transition nl n'l ', which we denote 

&nl—yn'l' • 

In our model atom, singlet and triplet terms are often 
separated, and thus we would like to estimate how this to¬ 
tal cross-section is distributed among final spin states. We 
considered a representation where the spin of the Rydberg 
electron is coupled to the spin of the other electrons in 
the target atom S c (assuming this to be a good quantum 
number) to give the total electronic spin of the Rydberg 
atom, 5. Recoupling of angular momenta gives the required 
spin-changing scattering amplitude in terms of the singlet 
and triplet scattering amplitudes (i.e. in the representa¬ 
tion in which the Rydberg electron spin and the hydrogen 
atom spin are coupled), and the differential cross-section 
for change of spin S —> S', S ^ S', is found to be 

(25" + 1) 

“ 2(25 c + l) exc?i [ ) 


where a exc h is the differential cross-section f or e l ectr on 
exchange in el astic e _ +H collisions (e.g. iFie 3 119581 
Burke fc Schevl 19621). T he derivation essentially follows 


Dalgarno fe Rudeel (|1965l ), except that in our case the elec¬ 


tron spin is coupled to the core electronic spin and not to 
the nuclear spin. We note that S' = 5 C ± 1/2 and the factor 
(25' + l)/2(25 c + 1) is always in the range 0 to 1 and rep¬ 
resents the probability that an electron exchange leads to a 
change in the total electronic spin of the Rydberg atom. In 
the limit p —1 0, —/ ± —» a , where a ± are the scattering 
lengths, which have values of a + = 5.965 and a~ = 1.769 
atomic units (iSchwartdll96 1). In this limit, for the case of 
S c = 1/2, via comparison of the expressions for the cross 
sections we find that 


^nljS—O—^n'l'jS'—O = 0.706 G n i^ n 'i> , 

(9) 

&nl,S=Q—tn , l',S'=l = 0.294 CT n /^. n ///, 

(10) 

^nljS—l—^n'l^S'—O — 0.098 

(11) 

&nl,S=l— >-n , / , ,iS' / =l — 0.902 dni—^n'i' . 

(12) 
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Mg i 



Mg ii 


3s-1S 
3p-3P* 
3p-1 P* 
4s-3S 
4s-1S 
3d-1D 
4p-3P* 
3d-3D 
4p-1 P* 
5S-3S 
5s-1S 
4d-1D 
4d-3D 
5p-3P* 
4f-1 F* 
4f-3F* 


20.9 ■ 

19.8120.9 
18.7119.8 
17.6118.7 
16.5117.6 
15.4 ■ 16.5 

14.3 M 15.4 

13.2114.3 
12.1 Hi 3.2 
11 .oll 2.1 

9.9 □ 11.0 
8.8 ■ 9.9 
7.7 ■ 8.8 
6.6 ■ 7.7 
5.5 EH 6.6 
4.4CH5.5 
3.3CH4.4 
2.2D3.3 

1.1 D 2.2 

□ l.l 


16.1 ■ 

15.3 ■ 16.1 
14.5H15.3 
13.6 ■ 14.5 
12.8B13.6 
11.9 ■ 12.8 
11.0B11.9 
10.2 El 11.0 
9.4 ■ 10.2 
8.5IB9.4 

7.608.5 

6.807.6 
5.9 IE 6.8 

5.1 IB 5.9 

4.215.1 
3.4D4.2 
2.6D3.4 
1.7D2.6 
0.9 □ 1.7 

□ 0.9 


Fig. 5: Collision matrices for Mg I and Mgll at logr 500nm = —0.23 in the solar (3 D 0 )id atmospheric model, which 
corresponds to a temperature of 6000 K. The legend in each case defines a range of log C\j , where C\j is the total collision 
rate in s -1 , such that darker squares represent higher rates. Note the 3 p 2 level of Mg I is not collisionally coupled to 


other levels; see Sect. 13.31 
3.4. Test model atoms 

To understand the uncertainties in the modelling, we 
constructed various model atoms in which sets of colli- 


sional data were changed or removed compared to the final 
model. For ease of discussion, we adopted the following 
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abbreviation^ for the various collision processes: 

CE Collisional (de) excitation by electrons 

Cl Collisional ionisation by electrons 

CH Collisional (de)excitation by hydrogen 

CHO Charge transfer with hydrogen. 

For reference, our final model adopted for CE data RM 
data wherever possible, IP data for all remaining allowed 
transitions, and we used the recipe described in Sect. 13.3.11 
and defined by Table [5] for forbidden transitions (hereafter 
referred to as the “Y” recipe). The final model, which we 
give the label F, also includes data for Cl, CH and CHO as 
described in Sect. 13.31 

The first additional test model represents a basic model 
without any of the new data introduced here. This model 
contains only CE and Cl, and no hydrogen collision pro¬ 
cesses. For CE it employs IP data for all allowed transitions 
and uses the T recipe for all forbidden transition^. Thus, 
we refer to it, where appropriate, as the basic model, model 
B. 

The models F and B therefore represent two extremes, 
and it is interesting to investigate models that lie between 
these cases. In particular, we constructed models B+RM 
(introduces RM), B+RM+CH (further introduces CH), 
B+RM+CHO (introduces RM and CHO), B+H (introduces 
CH and CHO), which may be used to explore the relative im¬ 
portance of RM, CH, CHO, or hydrogen collisions in total, 
H, in causing the differences between model B and model 
F. To isolate the effects of spin changing collisions, we con¬ 
structed versions of the final model where spin changing 
collisions (i.e. between singlet and triplet states in Mgl) 
are removed for collisions with electrons (F—SCE), hydro¬ 
gen (F—SCH) or both (F—S). Transitions between levels of 
different spin can only be removed when the spin of each 
level is defined, but not in merged super levels. 

Finally, we constructed slightly modified versions of 
model F, in which the IP data were replaced with vR data 
and/or the Y recipe was replaced by a typical approxima¬ 
tion of setting the collision strength to unity, fl = 1. These 
models are labelled F(vR,T), F(vR,H = 1), F(IP,D = 1), 
noting that F(IP,T) would be equivalent to the final model 
F in this notation, and these models can be used to inves¬ 
tigate the impact of various approximations for CE among 
high-lying levels not covered by RM data. 

Table [7] gives an overwiew of these models. By examin¬ 
ing all of these models, we can determine for a given model 
atmosphere which of the newly introduced data or processes 
has produced the most important changes in each line, and 
this is reported in the last four columns of Tabic [T] and 
is part of our discussions of below, where we discuss the 
results for the various groups of lines. 


4. Comparisons with observations 

We now test our modelling by comparison with observed 
spectra. First, we compare our results with optical spectra 
of standard benchmark late-type stars, including the Sun. 

6 These notations are adopted from MULTI. 

7 The Y recipe is a new addition, based on new RM data. How¬ 
ever, its importance is limited and we chose to use it in the basic 
model to limit the number of degrees of freedom in our tests. Its 
effects are tested via the F(IP,fI = 1) model described later in 
this section. 


Table 7: Description of the collisional data for the different 
versions of model atoms described in Sect. Eli 


Name 

CE 

CH 

CHO 

Allowed 

For bidden 

B 

ip 

T 



B+RM 

RM+IP 

RM+T 



B+RM+CH 

RM+IP 

RM+T 

/ 


B+RM+CHO 

RM+IP 

RM+T 


/ 

B+H 

IP 

T 

/ 

/ 

F = F(IP,T) 

RM+IP 

RM+T 

/ 

/ 


F-SCE 

RM+IP 

RM+T—spin 
exchange 

/ 

/ 

F-SCH 

RM+IP 

RM+T 

No spin 
exchange 

/ 

F-S 

RM+IP 

RM+T—spin 
exchange 

No spin 
exchange 

/ 


F(IP,f!=l) 

RM+IP 

RM+(H = 1) 

/ 

/ 

F(vR,T) 

RM+vR 

RM+T 

/ 

/ 

F(vR,r2=l) 

RM+vR 

RM+(H = 1) 

/ 

/ 


Notes. RM+x means RM when available, otherwise method x. 


Then we compare them with the IR emission lines in the 
solar spectrum. 


4.1. Optical absorption lines 


We tested the performance of our non-LTE modelling using 
high-quality spectra of six stars with reasonably well-known 
fundamental parameters. This sample of benchmark stars 
was designed to cover a wide parameter range and includes 
the Sun, Procyon, Arcturus, and three very metal-poor 
stars ([Fe/H] = —2.5... — 2.0); a turn-off star (HD 84937), 
a sub-giant (HD 142083), and a red giant (HD 122563). For 
the Sun and Arcturus, spectra from the Kitt Peak ob¬ 
serva tory with R & 150 000 and S/N « 1000 are avail¬ 
able dHinkle et al.l 12000 b For the other stars, we retrieved 
VLT/UVES spectra at R « 47 000 from the POP archive 
of bright field stars ( Bagnulo et a lJ l2003h . 


Accurate effective temperatures for Procyon, Arcturus, 
and HD 122563 were derived from direct measurements of 
their ang ular d iam eters (see adopted parameters in Table 
l51). ICreevev et all (|2014f) presented fundamental parameters 
of HD 140283, although T e g is uncertain, and for HD 84937 
no such measurements have been publishe d, so we r elied 
inste ad on the infra-red flux method (Casagrandc et al. 
l2010 h For all stars except for the Sun, HiPPARCOS paral- 


itv (Ivan Leeuwen 

2007) 

from Jofre et al. ( 

20 + 


). Reference metalliciti es wer e taken 

_ _ __ ■■■_. _ for Arcturus and iBergemann et al.l 

(|2012h for all other stars. Microturbulence values were also 
taken from the literature studies and held fixed in the line 
formation calculations. 


We obtained ID and (3D) id model at mospheres fo r 
each star by interpolati ng in the MAR.CS (iGustafsson et al.l 
I2008H and STAGGER (iMagic et al.ll2013h grids. Minor ex¬ 
trapolation in the (3D) id grid was sometimes necessary to 
produce the desired model parameters. The I D spectrum 
synthesis code SME (IValenti fe Piskunov! 119961 . and subse¬ 
quent updates) was used to compute synthetic spectra in¬ 
cluding blending lines. To synthesise the lines in non-LTE, 
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8805 8806 8807 8808 

Wavelength [A] 


8805 8806 8807 8808 

Wavelength [A] 


8805 8806 8807 8808 

Wavelength [A] 


Fig. 6: Observed spectra of the non-LTE sensitive Mg I 8806 A line (pluses). The best-fitting (3D)idLTE (red dashed 
line) and non-LTE profile (blue dotted line) are overplotted and the corresponding Mg abundances stated in each panel. 

Table 8: Stellar parameters and average Mg abundances, including line-to-line la dispersions, for the benchmark stars. 


Name 

T e S 

[K1 

log g 

[Fe/H] 

^mic 

[km/s] 


ID 



(3D) 

ID 


LTE 

non- 

-LTE 

LTE 

non- 

LTE 

Sun (“> 

5777 

4.44 

0.00 

1.09 

7.57 

0.08 

7.57 

0.08 

7.67 

0.07 

7.66 

0.07 

Procyon ^ 

6545 

3.99 

-0.03 

2.00 

7.47 

0.04 

7.45 

0.04 

7.60 

0.05 

7.56 

0.06 

Arcturus( b ) 

4247 

1.59 

-0.52 

1.63 

7.42 

0.07 

7.33 

0.07 

7.52 

0.10 

7.38 

0.06 

HD 84937 ^ 

6408 

4.13 

-2.03 

1.40 

5.73 

0.04 

5.75 

0.03 

5.78 

0.05 

5.80 

0.03 

HD 140283 ( “) 

5777 

3.67 

-2.40 

1.23 

5.36 

0.10 

5.39 

0.10 

5.35 

0.16 

5.45 

0.09 

HD 122563 (b ) 

4608 

1.61 

-2.64 

1.50 

5.29 

0.13 

5.26 

0.13 

5.28 

0.08 

5.27 

0.08 


Model parameters taken from _Ber 2 emann_et_aL (2012), ^ Jbfre_e^_aL (2014) 


grids of pre-computed departure coefficient^ bi from MULTI 
(using model F) were read in and interpolated before they 
were applied to the LTE level populations computed by SME. 
For all interpolation and line-formation, the continuum op¬ 
tical depth at 5000 A was used as reference scale, which 
is consistent with the reference scale used when computing 
the spatial and temporal average of the full 3D models. We 
compared the spectral line profiles calculated by MULTI 
and SME used in this fashion, and they give practically 
identical results. 


8 Defined a s bi = m/n* where m is the non-LTE population 
and n* the LTE population of level i. 


Between eight and eleven Mg lines, selected to be rela¬ 
tively blend-free and situated between 4000 — 9000A were 
synthesised per star. Whenever possible, the best-fitting 
abundance was found by simultaneously varying the Mg 
abundance and the intrinsic line broadening, such as caused 
by macro-turbulence and rotation. It was often found that 
neglecting non-LTE effects in the cores of strong lines pre¬ 
vented a match between the LTE synthetic spectrum and 
the observations without artificially reducing all broaden¬ 
ing, even the instrumental, to unphysical values. In such 
cases, we instead fixed the broadening parameters to the 
values found in non-LTE and disregarded the line core when 
optimising the fit. By doing so, the agreement between LTE 
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and non-LTE abundances is typically better than if the 
same equivalent widths had been enforced, for instance. For 
lines weaker than 150 mA, the LTE and non-LTE line pro¬ 
files are similar, and a good fit to the full observed line 
profile can always be found. Sample (3D) id fits are seen in 
Fig. [6] for Mg I 8806 A. 

The choice of model atmosphere in ID or (3D) id and 
line-formation method in LTE or non-LTE affects the de¬ 
rived Mg abundances to various degrees and in various di¬ 
rections. Average abundances and line-by-line ler disper¬ 
sions are presented in Tabled]. We stress that our main goal 
is to illustrate the impact of the modelling method on a va¬ 
riety of different lines, not only those that are most suitable 
for abundance analysis. The average abundances we present 
should therefore not be taken too literally. For the Sun, non- 
LTE effects are very weak, of the order of 0.01 dex, while 
(3D)id abundances are higher by 0.05 — 0.15 dex than the 
corresponding ID abundances because of the hotter tem¬ 
peratures (~ 1%) in the line-formation region. Only two 
lines overlap with the recent study bv IScott et al. (201 l i. 
but the average solar abundances nevertheless match well 
(A(Mg) Scottet . al . =7.59±0.04). 

In Arcturus and Procyon the effects work in the same 
direction as in the Sun, but the sensitivity to non-LTE is 
higher. The over-population of Mg I penetrates to deeper 
layers and abundance corrections reach up to —0.1 dex in 
Procyon and up to —0.3 dex in Arcturus. In the metal- 
poor stars, most lines are subjected to a similar weaken¬ 
ing, that is, abundance strengthening, in (3 D)id LTE as 
in ID LTE because of the somewhat hotter temperatures 
at T 500nm « 1. However, lines that form relatively far out, 
such as the intercombination line 4571 A and the near 
IR 8806 A, are affected by the strong cooling of the sur¬ 
face layers that is characteris tic o f meta l-poor rad iatio n- 
hydrodynamical models (lAsplund fc Garcia Perezl 1200 11. 
Consequently, (3D) id LTE abundances are significantly 
lower for these lines, up to —0.25 dex in HD122563. How¬ 
ever, the effects of the surface cooling on line strengths 
is counteracted by non-LTE effects in the opposite direc¬ 
tion. The behaviour is quali tativel y similar to that found 
for resonance line s of C al in lLind et al.1 (l2013h and Fel in 
iBergemann et all (|2012f) . As also described in these stud¬ 
ies, the steeper temperature gradients of metal-poor (3D) id 
models than in ID models strengthens the overionisation of 
neutral species, such as Mg I, resulting in positive non-LTE 
corrections for all lines. Reassuringly, the line-to-line abun¬ 
dance scatter in (3D) id non-LTE is overall an improvement 
with respect to ID LTE. 

In Fig. |7] we illustrate how abundances of the bench¬ 
mark stars compare when limited to the diagnostics avail¬ 
able to the GALAH and Gaia-ESO survey for metal-rich 
and metal-poor stars. Evidently, systematic uncertainties 
are comparable to or even exceed the targeted precision (0.1 
and 0.05 dex for Gaia-ESO and GALAH, respectively). The 
subgiant and dwarf stars in our sample do not appear to 
suffer from strong differential effects, and sufficiently precise 
relative abundances can probably be obtained also with tra¬ 
ditional modelling of Mg, although they may not necessar¬ 
ily be accurate in an absolute sense. However, giants show 
strong differential non-LTE effects (0.1 — 0.2 dex), which 


9 The abundance has the usual definition A (Mg) = 
log {NmJNh) + 12, where Wig is the number density of mag¬ 
nesium atoms and ions. 
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Fig. 7: Top panel: Mg abundances of stars in the solar 
neighbourhood affected by different types of modelling. The 
spectral lines used are those covered by the GALAH sur¬ 
vey: 571lA and 769lA. The bottom panel shows the same 
for metal-poor halo stars, based on lines observed in the 
Gaia-ESO survey: 5528A and 8806A. 

must be corrected for to allow comparisons to less evolved 
stars. 

4.2. Infrared emission lines 

As discussed in Sect. [21 the formation of IR Mg I emission 
lines in the solar spectrum, discovered in the early 80s, 
took a decade and considerable effort to understand. Two 
main hypotheses arose: they form in the chromosphere un¬ 
der LTE conditions, or they form in the photosphere with a 
non-LTE source function increasing outwards. This debate 
was finally and unambiguously resolved bv iGarlsson et all 
(H993) in favour of a non-LTE line formation scenario of 
photospheric origin. Their modelling made several reason¬ 
able, yet ad ho c , assumptions rega r ding c ollision rates. As 
pointed out bv lRutten k. Carlssonl (I1994D . the 12 pm lines 
have non-LTE emissions peaks because of the strong colli- 
sional coupling between the pertinent levels. Based on the 
expectation that hydrog e n collisions wou ld provide such 
high rates (lOmontJ I1977T) . IGarlsson et all introduced high 
electron collisional rates for Z-changing collisions among Ry¬ 
dberg states to ensure detailed balancing between levels 
of equal n. A t temp ts at modelling since that t ime, such 
as IZhao et al.l (Il998l) and ISnndavist et aD (1 20081 ). have de¬ 
scribed hydrogen collisions using the highly debated Drawin 
formula scaled by a free parameter chosen to best fit ob¬ 
servations. Importantly, however, these studies showed that 
the IR. lines in the Sun and in red giants such as Arcturus 
are rather sensitive to the description of hydrogen colli¬ 
sions, and as these lines probe the highly excited parts of 
the atom, they are important tests of our ab initio mod¬ 
elling. 

Figure [8] compares theoretical spectra for the 7.3, 12.2, 
and 12.3 pm emission lines with the observed solar spec¬ 
trum across the solar disk. The synthetic spectral lines are 
calculated for various model atoms, using the curren tly 
accep ted solar photospheric Mg abundance (IScott et al.l 
l2014j) . A(Mg)=7.59, and employing the (3 D)id solar model 
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atmosphere. We note that our modelling does not include 
pressure-induced line shifts, and other atomic data for such 
shifts do not exist. Thus, we have no expectation to be able 
to model the asymmetries observed in the wings of these 
lines. 

First, we note the well-known and practically obvious re¬ 
sult that the lines cannot be reproduced in a photospheric 
model using LTE line formation. Second, we see that our 
model F reproduces the observations both in the core and 
the wings rather well, but certainly not perfectly. However, 
it certainly reproduces the observations better than other 
model atoms of interest. As noted in Sect. 13.3+1 on physical 
grounds the vR formula is expected to overestimate the CE 
rates because it is based on the Born approximation. The 
vR formula results in rates higher than those given by the 
IP method, producing stronger collisional coupling between 
levels of different n and thus reducing the strength of the 
emission. The predicted line cores for the model F(vR,T) 
are in all cases far too weak. In addition, we also see that 
a model atom with no hydrogen collision processes, model 
B+RM, produces line profiles that are generally too strong, 
both in the core and in the wings. This supports the conclu¬ 
sion of earlier studies that additional collisional mechanisms 
are required. Thus, we see that our model F employing the 
IP method for CE rates among Rydberg states, and the 
Kaulakys method for CH among Rydberg states, generally 
produces the line profiles better than models using the vR. 
method or ignoring CH. The fact that the astrophysical 
comparison agrees with expectations based on the physical 
characteristics of the vR and IP methods, gives us confi¬ 
dence that our ab initio modelling of collision rates among 
Rydberg states due to CE and CH is satisfactory. 

We note that the emission peaks are sensitive to macro¬ 
turbulence (see Fig. [8]) and the wings can be better fitted 
with a change in abundance of order 0.1 dex. As might be 
reasonably expected, the CHO processes do not significantly 
affect these lines because CHO processes only involve rela- 
tiv ely low-lying stat e s, and thus any effect must be indirect. 

ISundavist et al. d2008ih studied these lines in a small 
sample of K giants including Arcturus, in which they 
showed that the 12.2/um line was particularly sensitive to 
the modelling of hydrogen collisions. This is due to the small 
influence of electron collisions in such atmospheres. Fig¬ 
ure El compares our results with their observations (kindly 
provided by Nils Ryde). The LTE calculation is again not 
able to form this line in emission. For our model atom F, 
the match in the wings is good, while the emission peak is 
stronger than in the observed spectra; however, we consider 
the agreement reasonable given, for example, the uncertain¬ 
ties in stellar parameters and that the macroturbulence and 
abundance have not been fine-tuned. As is the case for the 
Sun, our model F reproduces the line better or as well as 
the other models of interest. In particular, neglect of hydro¬ 
gen collisions, model B+RM, again produces a peak con¬ 
siderably too strong. We note that in this case, the F(vR) 
model performs marginally better than the F(IP) model; 
however, the differences are small and we judge that this is 
insignificant with respect to the other modelling uncertain¬ 
ties mentioned above. 

5. Spectral line behaviour 

In this section we examine the behaviour of the modelled 
spectral lines, in particular the influence of different colli- 
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Fig. 9: Profile of the 12.2/im line observed in Arcturus by 
Sundqvist et al. (black circles) and compared with our cal¬ 
culations using the (3D) id atmospheric model and the non- 
LTE abundance taken from Table [5] The black lines show 
synthetic profiles using different model atoms (see Table Q, 
and are all convolved with a macroturbulent velocity of 5.2 
km/s and a rotational broadening of v sin* = 2.1 krn/s. The 
grey line is the profile using the F model atom but without 
macroturbulence. 

Table 9: Stellar parameters of the atmospheric models used 
in Sec [5] 


Name 

T e fi 

log g 

[Fe/H] 

^mic 


[K] 

[cm/s 2 ] 

[dex] 

[km/s| 

Dwarf Poor 

6000 

'4.5 

-2.0 

1.0 

Dwarf Rich 

6000 

4.5 

0.0 

1.0 

Giant Poor 

4500 

1.0 

-2.0 

2.0 

Giant Rich 

4500 

1.0 

0.0 

2.0 


sional processes on line strengths and the abundances that 
would be derived from the lines. Tests were performed on a 
small grid of four ID model atmospheres from the MARCS 
grid <|Gustafsson et all I2008T) with parameters described in 
Table [2] where u m i c is the micro turbulent velocity used 
in calculating the model atmospheres. The same value is 
used in the line formation calculations. Additionally, we 
performed a number of test calculations for the solar case, 
using the MARCS solar model. We also performed calcula¬ 
tions using the photos pheric r eferen ce model with a chro¬ 
mosphere from iMaltbv et all (|1986f) . hereafter MACKKL, 
and the ID-averaged 3D model for the Sun, (3Dq)id, from 
iMagic et all (|2013f ). Furthermore, to investigate possible ef¬ 
fects of the chromosphere, we also made a hybrid model 
using MACKKL at rsoonm > — 3, while replacing the upper 
layers with those from the (3Dq)id model. Where com¬ 
parisons are made with solar spectra, we used a Gaussian 
macro-turbulence with u mac = 2.0 km/s. 

As we are most interested in effects on the derived abun¬ 
dance of Mg, we studied the influence of the various colli- 
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Fig. 8: Solar IR line profiles at various positions on the disk. The first row is the 7.3 /xm line, the second row the 
12.2 /xm line, and the fi nal row the 1 2.3 /xm line. The columns are /x = 1.0, 0.5, and 0.2, in t hat o rder. The circles are 
observations taken from I Chang et al.l (|l99lh for the 7.3 /xm line and from iBrault fe Novesl (|l983h for the two 12 /x m 
lines. The synthetic spectral lines are calculated for various model atoms using A(Mg)=7.6 and employing the (3D)i£> 
solar model atmosphere. For the synthetic profiles, a macroturbulent velocity of 3.15 km/s and a rotational broadening of 
v sinx = 1.6 km/s are used. For model F, the profile without macroturbulence is also shown (grey full line) to demonstrate 
its effect. 


sion processes on the line formation of Mg lines through the non-LTE corrections, defined by 

A^(Mg) NLTE _ LTE = A(Mg) NLTE - A(Mg) LTE . 


(13) 
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The abundance corrections discussed in this section are 
based on matching equivalent widths. We note that they 
therefore cannot be compared to the changes in abundances 
seen in the last section, where profile fitting was used. 

For the lines to be discussed here, various non-LTE 
mechanisms interact to varying degrees to produce the final 
line profiles. Generally, overionization leads to an underpop¬ 
ulation of the low-lying levels of Mg I, and a number of Mg 
lines lead to photon losses, resulting in depart ures f rom LTE 
also among intermediate levels (e.g. ICarlsson et aljfl992l h 
These radiative mechanisms may be counteracted or even 
propagated by collisional processes. Charge exchange pro¬ 
cesses (CHO) drive the populations of specific intermediate 
levels, in particular those around 4s 1 S, toward detailed bal¬ 
ance with the ground state of Mg II, which is generally close 
to LTE populations. CH couples low-lying and intermedi¬ 
ate levels propagating departures from low-lying to higher 
levels. 

5.1. Mg I lines 

3829, 3832, and 3838 A (3d 3 D-3p 3 P°): The UV triplet 
is one of the strongest features of Mgi and is often used 
for abundance determinations in metal-poor stars. In solar 
metallicity stars these lines are significantly blended, and 
also sensitive to chromospheric effects. The abundance cor¬ 
rections and their sensitivity to collisional processes are pre¬ 
sented in Fig. HOI Model B leads to very large corrections for 
these lines. However, in model F abundance corrections are 
general quite small. In dwarfs the corrections are practically 
non-existent, but if collisions with H were not included, we 
would find corrections of ~ 0.2 dex in the metal-poor case. 
Thus, in dwarfs CH and CHO play an important role in 
the abundance correction calculation. In giants, corrections 
are also small, 0.05 dex for abundance values of inter¬ 
est, and neglecting CH leads to much larger corrections, 
of about 0.2 dex in the metal-poor case. This means that 
in giants, CH stronglycontributes to the abundance correc¬ 
tion and CHO has a negligible effect. If hydrogen collisions 
are included, the use of more approximate models for CE, 
such as F(vR, fl = 1) or F(IP, fi = 1) does not change 
the results significantly. Spin change transitions have only 
a weak effect in the two dwarf stellar models, the correction 
changing by ~ 0.05 dex in the metal-poor case for the F—S 
model. In giants, the F—S model changes the corrections by 
~ 0.07 dex with the relative effect of SCE being strongest 
at solar metallicity and the effect of SCH strongest at low 
metallicity. 

4167, 4703, 5528, and 8806 A (nd 3 D-3p 1 P°): For a given 
stellar model, all four of these lines behave similarly, but the 
effects becomes stronger at longer wavelengths, where the 
lines form at shallower depths with larger population de¬ 
partures. Figure [TU] shows the abundance corrections for the 
5528 A line. Model B leads to large corrections in dwarfs, of 
about 0.2 dex at the abundances of most interest. The cor¬ 
rections with model F are, however, small. The results for 
4167, 4703, and 5528 were most sensitive to the introduc¬ 
tion of hydrogen collisions generally (both CH and CHO), 
while for the 8806 line CH was particularly important. This 
is because the populations of the levels involved are simi¬ 
larly affected by the introduction of CHO; the upper level 
of the 8806 A transition, 3d 3 D, is low-lying and both lev¬ 


els are thermalized by CHO to similar degrees in the line 
formation region, thus leaving S„ unaltered. 

4571 A (3p 3 P°-3s 1 S): The intercombination line, de¬ 
spite its weak oscillator strength, is observable across a 
wide range of stellar parameters because it involves the 
ground state of Mgi. The abundance corrections are shown 
in Fig. [TOJ In dwarfs, at abundances of typical interest, this 
line is formed close to LTE conditions as the collisional cou¬ 
pling of the two levels in different spin systems dominates 
the weak radiative coupling. In giant models where colli¬ 
sional coupling is significantly decreased, this line shows 
significant departures from LTE, of the order of 0.2 dex at 
solar metallicity and almost 0.4 dex in the metal-poor gi¬ 
ant model. In metal-rich giants, CHO is much weaker than 
either CH and RM. This line is expected to be sensitive to 
spin change transitions. CE rates between the upper and 
lower levels are 10 8 times higher than the corresponding 
CH rates, making this line most sensitive to spin change 
CE rates even in very metal-poor stars. 

5167, 5172, and 5184 A (4s 3 S-3p 3 P°): The abundance cor¬ 
rections for one of the Mgi b triplet lines are shown in 
Fig. [lOl Model B gives large corrections for these lines, but 
model F gives corrections that are very small for dwarfs, 
and quite small for giants, lower than 0.05 dex. As with 
the 8806 A line, the departure coefficients of the upper and 
lower levels are similarly affected by CHO, and thus CHO 
has little effect on the line strength. The main change is due 
to the introduction of CH. CH decreases the population of 
the upper level in the line formation region and increases 
the populations of the lower levels, bringing them closer to 
relative LTE. 

8736 A (7/ 3 F°-3d 3 D) and 8710, 8713, and 8718 A 
(7d 3 D-4p 3 P°): These lines are n = 7 -o- n' = 3 tran¬ 

sitions involving triplet states, that is, between low-lying 
and high-lying levels. The abundance corrections for one 
of the lines are shown in Fig. [101 These lines show exam¬ 
ples of behaviour where the addition of collision processes 
can drive away from LTE, rather than towards it. As these 
levels lie approximately 3-4 eV from the continuum, they 
are subject to overionisation f r om u ltraviolet radiation. In 
the solar case, ICarlsson et al.1 (119921) have shown that the 
triplet 3d and 4 p levels are among those involved in transi¬ 
tions that become optically thin in the photosphere, which 
in turn leads to photon losses and thus photon suction. 
Now, together with these radiatively driven non-LTE pro¬ 
cesses, the CH and CHO processes involving these states are 
among the strongest. 

The complex interplay of these processes leads to the 
situation seen in Fig. 1101 where rather than the more usual 
convergence from model B towards the F model, the in¬ 
troduction of additional collisional processes has less pre¬ 
dictable effects. For example, in the metal-poor dwarf case, 
with model B, rather large corrections are seen. Introduc¬ 
tion of the RM data via model B+RM leads to slightly 
increased corrections. Then, contrary to the naive expec¬ 
tation that additional collisions drive towards LTE, the in¬ 
troduction of CH, model B+RM+CH actually increases the 
abundance corrections. The left panel of Fig. fill shows de¬ 
parture coefficients of the levels involved in the formation of 


Article number, page 18 of 1221 








Y. Osorio et al.: Mg line formation in late-type stellar atmospheres: I. The model atom 


A(M S ) LTE 

383.8 nm 

5.5 6.0 6.5 7.0 7.5 8.0 



A(Mg) LTE 

552.8 nm 



5.5 6.0 6.5 7.0 7.5 8.0 


A(Mg) LTE 


A(Mg) LTE 


457.1 nm 
5.5 6.0 


A(Mg) LTE 
6.5 7.0 



518.4 nm 
5.5 6.0 


A(Mg) LTE 
6.5 7.0 





0.4 




0.3 








0.2 



AA(Mg) 

0.1 




0.0 

0.5 



-0.1 

0.4 

g 



0.3 




0.2 

AA(Mg) 



0.1 




0.0 





—-__ 

F- 

B+RM+CH. 

B+RM+CH0- 

B+RM- 

B- 




_..- "A 




-= 

. _ — 


...... 






6000/4.5/0.0 



6000/4.5/-2 

.0 


— 


- - 


4500/1.0/0.0 



4500/1.0/-2 

.0 


A(Mg) LT 


6.5 7.0 

A(Mg) LTE 


A(Mg) LTE 

871.8 nm 

5.5 6.0 6.5 7.0 7.5 8.0 



A(Mg) LTC 


Fig. 10: Abundance corrections for the 3838 A, 5528 A, 4571 A, 5184 A, and 8718 A lines in four different atmospheric 
models and using modified model atoms as described in the text. The model parameters are given at the bottom left of 
each figure with T e ff/log ^/[Fe/H]. The shaded region indicates the likely abundance of Mg for a normal star of given 
metallicity. The black circles show the calculated abundances; some lines have fewer points because the calculations with 
those parameters did not converge. Note that the model atoms shown in the 4571 A line case differ from the others. 


the 8718 A line in a metal-poor giant atmosphere, where we 
see that the introduction of CH leads to populations of the 
4 p level at around log(r 5 oonm) = —2 to —3 closer to the LTE 
population, but the population at around —1, where the line 
forms, is in fact farther from LTE. This occurs because deep 


in the atmosphere, CH processes lead to increased coupling 
with the lower-lying levels, especially 4s and 3 p levels, which 
are subject to strong photoionisation and thus violent ove¬ 
rionisation, leading to underpopulation. However, we see in 
this model atmosphere that the dominant change is due to 
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T 
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Fig. 11: Departure coefficients b = n/n LTE of the levels involved in the 8718 A line, 7d (blue) and 4p (red). The left 
panel shows the metal-poor dwarf case, the right panel the metal-poor giant case. The insets show the corresponding line 
profiles. Note that in the giant case, abundance corrections corresponding to some models for which b are shown here 
are not shown in Fig. [TO] because the models did not converge for at least three abundance values. 


the introduction of CHO, which increases the coupling be¬ 
tween the 4s, 3 d and 4 p levels with the continuum and thus 
brings the population closer to relative equilibrium with Mg 
II and thus towards LTE populations, particularly in the 
line-forming region. For the metal-poor giant model, with 
departure coefficients shown in the right panel of Fig. fTTl 
model B shows an overpopulation of the 4 p level in the line¬ 
forming region log(r 5 oo n m) = —1 caused by photon losses. 
The introduction of CH processes again pulls the popu¬ 
lations towards those of the lower-lying levels, which are 
underpopulated between log(r 5 oo n m) = —1 and 0.5. This 
causes an underpopulation at large depths, but brings the 
population to near LTE levels in the line-forming region. 
Introduction of CHO processes brings the population deep 
in the atmosphere back towards LTE, which in turn rein¬ 
troduces a slight overpopulation in the line-forming region. 

5.2. Mg II lines 

Since Mg II is the dominant ion of Mg in late-type stars, 
non-LTE effects are expected to be weaker in Mg II than in 
Mg I lines. In the four test model atmospheres studied here, 
the abundance corrections all had magnitudes lower than 
0.1 dex for the lines we analysed (see table [1]). 

6. Conclusions 

We have presented a new model atom for Mg suitable for 
non-LTE line formation studies in late-type stars. It in¬ 
cludes data for electronic collision rates calculated in this 
paper by the i?-matrix method, and we have further used 
these results to formulate a method to estimate transition 
rates for forbidden transitions not covered by the i?-matrix 
calculations, and not possible to be estimated by the vR 
and IP methods. Recent data for hydrogen collisions, in¬ 
cluding charge transfer processes, calculated by some of us, 
were also applied. Some new data for collisional broaden¬ 
ing were calculated as well. The new model atom, used to¬ 


gether with ID model atmospheres in the context of stan¬ 
dard non-LTE modelling, compares well with observations. 
The modelled spectra agree well with observed spectra from 
benchmark stars, showing much better agreement with line 
profile shapes than LTE modelling. The line-to-line scatter 
in the derived abundances shows some improvements com¬ 
pared to LTE (where the cores of strong lines must often 
be ignored), particularly when coupled with averaged 3D 
models. The observed Mg emission features at 7 and 12 pm 
in the spectra of the Sun and Arcturus, which are sensitive 
to the collision data, are reasonably well reproduced. 

The spectral line behaviour and uncertainties were ex¬ 
plored by extensive experiments in which sets of collisional 
data were changed or removed on a small grid of theoret¬ 
ical model atmospheres. Charge transfer in collisions with 
H is generally important as a thermalising mechanism in 
dwarfs, less so in giants. Excitation due to collisions with H 
was found to be quite important in both giants and dwarfs. 
The i?-matrix calculations for electron collisions also lead 
to significant differences compared to when approximate 
formulas are employed. The modelling predicts non-LTE 
abundance corrections AA(Mg) NLTE _ LTE in dwarfs, both so¬ 
lar metallicity and metal-poor, to be very small (of about 
0.01 dex), even smaller than found in previous studies. In gi¬ 
ants, corrections vary greatly between lines, but can be as 
large as 0.4 dex. Our results emphasise the need for accurate 
data for electron and hydrogen collisions for precise non- 
LTE predictions of stellar spectra, but demonstrate that 
such data can be calculated and that ab initio non-LTE 
modelling without resort to free parameters is possible. 

Only two other elements, Li and Na, have collisional 
data for non-LTE modelling of comparable quality to that 
which has been used in this work. RM calculations have 
been performed for excitation of lo w-lyi ng levels of Li 
and Na by elect ron collisions by lOsorio et al.1 <|201lh and 
iGao et all (120101 ) resp e ctivel y. Before calculating the RM 
data, the data of iParkl (119711) were used in non-LTE mod¬ 
elling, and these data typically agree to within 80% of the 
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RM data, significantly better than the vR or IP methods. 
Thus, the new RM data for electron collisions did not have 
a strong effect on non-LTE modellin g, sinc e the existing 
data were r elativ ely good; see lOsorio et alJ (12011!) for Li, 
iLind et all (|2011h for Na. Data for hydrogen collision pro¬ 
cesses . including cha r ge tra nsfer, have been calculated for 
Li bvlBarklem et, al.l (|200.°I) . and for N a bvlB arklem et al.l 
( 2010T) . and used in non-LTE studies bv lBarklem et ahland 
ILind et all (1200911 in the case of Li and ILind et al. ( 201 Ilf 
for Na. These studies found for both elements that CH pro¬ 
cesses have practically no effect on the non-LTE modelling, 
while CHO, in particular that involving the first excited S- 
state, had a significant effect because it brought this level 
into closer detailed balance with the continuum. These are, 
however, alkalis with only one spin system. In contrast, Mg 
has two spin systems, leading to the added feature of inter¬ 
system transitions. Furthermore, many more lines of Mg 
are observable in stellar spectra from IR to UV. Our mod¬ 
elling has shown that the situation for slightly more com¬ 
plex atoms such as Mg is far more complicated than for 
alkalis, with the new data for CE, CH, and CHO all affect¬ 
ing at least some lines. 

We plan to use the current model atom to produce 
grids of departure coefficients and abundance corrections 
in ID and (3 D)id theoretical model atm ospheres, specifi¬ 
cally the MA RCS dGustafsson et al1l2008ll and STAGGER 
(jMagic et al.l 1201 .‘ill grids. Grids of abundance corrections 
can be directly applied to studies performed in LTE. How¬ 
ever, this depends on the use of equivalent widths, which is 
only strictly useful for weak lines. Grids of departure coef¬ 
ficients allow non-LTE corrections to be made at the spec¬ 
trum synthesis level, and thus is far more robust when using 
stronger lines. We have demonstrated that non-LTE calcu¬ 
lations are necessary for Galactic archaeology studies with 
large spectroscopic surveys of the Milky-Way. In particular, 
Mg abundances of giant stars suffer from strong differential 
non-LTE effects with respect to dwarfs and subgiants, com¬ 
parable to or exceeding the targeted precision. Using grids 
of departure coefficients computed with our model atom, 
on-going and future surveys can properly model these ef¬ 
fects. 

Finally, we note that a comparison with centre-to-limb 
variations of optical lines in the Sun would provide a 
strong complement to the tests of the modelling performed 
here. Such an observational campaign is underway with the 
Swedish Solar Telescope and will be the subject of future 
work. Ideally, this should be done in the context of full 3D 
non-LTE modelling. 
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